Bipolar charge-carrier injection in semiconductor /insulator /conductor 
heterostructures: self-consistent consideration 
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A self-consistent model of bipolar charge-carrier injection and transport processes in a semi- 
conductor/insulator/conductor system is developed which incorporates space-charge effects in the 
description of the injection process. The amount of charge-carriers injected is strongly determined 
by the energy barrier emerging at the contact, but at the same time the electrostatic potential 
generated by the injected charge-carriers modifies the height of this injection barrier itself. In our 
model, self-consistency is obtained by assuming continuity of the electric displacement and of the 
electrochemical potential all over the system. The constituents of the system are properly taken 
into account by means of their respective density of state distributions. The consequences resulting 
from our model are discussed on the basis of an indium tin oxide/organic semiconductor/conductor 
structure. The distributions of the charge carriers and the electric field through the electrodes and 
the organic layer are calculated. The recombination- and current-voltage characteristics are ana- 
lyzed for different heights of injection barriers and varying values of the recombination rate and 
compared with the measured current-voltage dependences for an indium tin oxide/poly(phenylene 
vinylene)/Ca structure. The voltage dependences of the recombination efficiency for the differ- 
ent values of injection barriers and recombination rate reveal optimum conditions for the device 
performance. 

PACS numbers: 73.40.Lq, 72.80.Le 
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I. INTRODUCTION 

Though the problem of bipolar charge injection in in- 
sulator media is many years o\dl££& } it anew became 
of great interest during last years in view of application 
of insulating materials as basic elements of electronic de- 
vices such as ferroelectric random access memories (FeR- 
AMs) 5 , organic field-effect transistors (OFETs) or or- 
ganic light-emitting diodes (OLEDs)£i&£. In the lat- 
ter devices organic semiconductors are used which show 
many properties of dielectric materials, especially rela- 
tively large band gaps and thus, a low intrinsic charge 
carrier density. The charge carriers have to be in- 
jected into the organic layer from the electrodes and 
thereby must overcome the injection barriers at the or- 
ganic/electrode interfaces. It is experimentally estab- 
lished that the injection conditions influence significantly 
the OLED performanc e) 10 ' 11 ' 12 ' 13 . 

Different models for the injection process are proposed 
in the literature. For low injection barriers, one expects 
the contact to be Ohmic, meaning that the contact is 
able to supply more charges per unit time than the bulk 
of the insulator can transport. In this case, a space- 
charge region is formed and the electric field at the in- 
terface is supposed to vanishii. Because excess charge- 
carriers dominate the charge transport in insulators, a 
space-charge-limited current (SCLC) density of the form 
j ~ V 2 /L 3 is observed (in the absence of charge-carrier 
traps), where L is the sample thickness and V is the 
applied voltage. Current-voltage (IV) characteristics of 
space-charge limited devices are determined by the bulk 
properties of the insulator with no influence of the con- 



tact propertie d 15 ' 16 ' 17 . 

For high injection barriers, one anticipates the injec- 
tion rate across the conductor/insulator interface to dom- 
inate the IV characteristic of the system. The models 
to describe injection are the Fowler-Nordheim (FN) tun- 
neling model or the Richardson-Schottky (RS) model for 
thermionic injection—. The FN model assumes tunneling 
through a triangular barrier into an unbound continuum 
of states. The RS model on the other hand describes 
charge injection as a thermally activated hopping over 
the potential barrier, where barrier lowering due to the 
superposition of the external electrostatic potential and 
the image-charge potential is considered. An attempt to 
incorporate space-charge effects into the hopping injec- 
tion model was undertaken in Refs. H9II20I . However, all 
these models consider injection as a single electron pro- 
cess. 

An alternative description of the injection process is 
given by the drift-diffusion theory involving electron- 
electron interaction in a mean-field approximation^!. 
In this homogeneous continuum model widely used for 
the description of conventional crystalline semiconduc- 
tors, the drift-diffusion equation in combination with the 
Poisson equation involves space-charge effects, but meets 
the problem of self-consistency in the boundary condi- 
tions. Namely, the electrostatic potential generated by 
the injected charge carriers modifies the injection barrier, 
but on the other hand the amount of charge-carriers in- 
jected per unit time depends on the barrier height. Thus, 
the values of the electrical field and charge-carrier density 
at the interface cannot be imposed but have to be found 
self-consistently implying the modification of the injec- 
tion barrier, too. An effort to perform a self-consistent 
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analysis of the contact phenomena was made in Ref. \2% 
However, the current and the field in this approach were 
not defined in a self-consistent way. This can be per- 
formed if the drift-diffusion theory involves both, the in- 
sulator and the conductor side of the system, in anal- 
ogy to a strongly asymmetric pn-junction. In the work 
of Neumann et al£^, a self-consistent numerical treat- 
ment of the unipolar injection and transport processes in 
a conductor/insulator/conductor device was presented in 
which continuity of the electrochemical potential and the 
electric displacement was assumed everywhere in the sys- 
tem, in particular at the contacts. The conductor and the 
insulator were characterized by their specific density of 
state (DOS) distributions. The exact analytical solution 
in a self-consistent manner for the charge injection across 
the sole conductor/insulator interface has been also ob- 
tained^. 

This model does not yet allow for a self-consistent de- 
scription of organic light emitting diodes, where bipo- 
lar transport takes place d 9 ' 25 ! 26 . To do so, one has to 
account for the recombination of charge carriers (holes 
and electrons) which are usually present in the organic 
constituent of the heterostructure under consideration. 
In such systems the recombination is often believed to 
be of Langevin typ o 27 i 28 , i.e. is considered as a bi- 
molecular process. Many OLED models involve the 
recombination of electrons and holes (see, for exam- 
ple, Refs. iQIISIlilMIlllllillllli). Probably, the most 
complete description was presented by Malliaras and 
Scot t 35 i 36 , who included a surface recombination at the 
contacts as well as a bimolecular recombination kinetics 
and diffusion in the bulk. However, the assumed bound- 
ary conditions for Ohmic contacts excluded diffusion at 
the contacts, where it is most important. A compre- 
hensive one-dimensional numerical model accounting for 
space-charge effect was developed by Tutis et al. which 
comprised the hopping transport in single and bilayer 
devices and tunnel injection from electrodes. However, 
besides the injection barriers given by the bare difference 
of the chemical potential in the metal and the frontier 
molecular orbital in the organic material this model in- 
volves a range of artifacts such as tunneling factors, effec- 
tive attempt frequencies, etc. Hence, in spite of a good 
agreement with experiments, this numerical tool is very 
complicated. 

In the present paper, we apply the diode model given in 
Ref. to elaborate a self-consistent OLED model. The 
properties of a semiconductor/organic/conductor struc- 
ture are modeled where holes and electrons are simulta- 
neously injected from an indium tin oxide (ITO) anode 
and a Ca cathode into the organic layer, respectively. 
The focus will be on the influence of the injection barriers 
and the recombination on the spatial distribution of the 
injected charge carriers and on the resulting IV charac- 
teristics as well as recombination efficiencies. Therefore, 
for simplicity, trap states and a field dependent mobility 
of the organic semiconductor constituent will be excluded 
at this stage. 



The paper is organized as follows. In Sec. II we present 
the theoretical model of the charge injection in semicon- 
ductor/insulator/conductor structure. Then, using spe- 
cific material parameters, we model the spatial distri- 
butions of the charge-carrier density and of the electric 
field (Sec. IIIA). In Sec. IIIB, the IV characteristics of the 
system (Sec. IIIB) are calculated varying the heights of 
injection barriers and the rate of the Langevin recombina- 
tion. We also discuss the calculated voltage dependences 
of the recombination efficiency of the considered struc- 
ture (Sec. IIIC). Our model is applied to a measured IV 
characteristic of an OLED composed of a poly(phenylene 
vinylene) (PPV) organic single layer sandwiched between 
ITO and Ca electrodes (Sec. HID). Finally, in Sec. IV our 
results are summarized. 



II. MODEL 

Let us consider an insulator of thickness L sandwiched 
in between a semiconductor and a conductor electrode. 
The insulator is supposed to be extended over the space 
with —L/2 < x < L/2, whereas the semiconductor and 
conductor electrodes are extended over the half-spaces 
with x < —L/2 and x > L/2, respectively. We assume 
that the semiconductor electrode injects holes into the 
insulator layer while the conductor electrode injects elec- 
trons. The band structure of the system under consider- 
ation is shown schematically in Fig. [TJ 

A. Electrodes 

Here we consider ITO as the hole-injecting electrode 
being an electron-conducting semiconductor with deep 
laying conduction band 21,38 . In the Thomas-Fermi ap- 
proximation, one can deduce the electrochemical poten- 
tial k s of the hole-injecting electrode as a function of the 
spatial coordinate x, 

K s (x) = ^— [3n 2 n s (x)] 2/3 - e<j) (x) , (1) 
2m* L J 

where n s (x) is the electron density in the electrode, m* 
is the effective mass of the charge carriers in the semicon- 
ductor, h is Planck's constant, e is the elementary charge, 
and (j> (x) is the electrostatic potential (notice that the 
energy level E = coincides with the bottom of the 
semiconductor conduction band). 

In general, everywhere in the heterostructure the elec- 
trochemical potential k (x) relates the steady-state cur- 
rent density j with the charge-carrier density. For a 
one-dimensional geometry, the current remains constant 
across the whole space and j is given by the conductivity 
a and the derivative of k (x) . 

Accordingly, for the hole-injecting electrode, 

/ \ dn s (x) . . 

j=fi s n s (x) — , (2) 
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where [i s is the charge-carrier mobility in the semicon- 
ducting electrode. The electric field F s (x) in the elec- 
trode obeys Gauss law, 



Sn s (x) 



(3) 



where e s is the relative permittivity of the electrode, 
8n s (x) = n s (x) 

~ n <x>.s is the excess electron density near 
the interface with respect to the electron density in the 
conduction band at infinite distance from the semicon- 
ductor/insulator interface, n^,^. The value of this excess 
density is supposed to be small in comparison with the 
background electron density, |<5n s (x)| <C n SiOC . Hence, 
the linearized Thomas-Fermi approximation can be ap- 
plied, leading to a differential equation for F s , 



^ = -l 2 TF , s F' s '(x)+F s (x), 



TF/j W T * s \*) , (4) 

with the conductivity in the electrode a s = e^ s noo,s and 



Itf,! 



' 2 £s e 0^oo,s 
3 C^^oc.s 



(5) 



being the Thomas-Fermi screening length of the elec- 
trode. Here, Koo iS is the Fermi level of electrons with 
respect to the bottom of the conduction band. 

Since the space-charge zone in the electrode is of finite 
thickness, the gradient of F s (x) has to vanish at infi- 
nite distance from the contact, and the solution for the 
electric field in the electrode reads 



F s (x) = i- 



, x + L/2 \ 
< x 'i' I i TFg ) ' ( 6 ) 



where the field in the electrode at the semiconduc- 
tor/insulator interface, F s (— L/2), is unknown and has 
to be determined by the boundary conditions. 

The distribution of the electric field in the electron- 
injecting electrode (conductor) is found in the same way 
as above. Even though this derivation was performed 
for a semiconductor it may be still meaningful for simple 
metals^. The electrochemical potential in the conductor 
reads 



boundary conditions, the conductivity a c — e^ c noc,c 
the Thomas-Fermi screening length 



TF,c 



2 Cc^0^oo,c 

3 C^77-oo,c 



and 



(9) 



Here \x c is the electron mobility of the metal electrode, 
Koo,c is the Fermi level of electrons with respect to the 
conduction band bottom, n.oo, c is the background elec- 
tron density in the conduction band at infinite distance 
from the conductor/insulator interface and e c is the rel- 
ative permittivity of the conductor. 



B. Insulator 

The energetic differences between the Fermi levels Kqo 
in the electrodes and the bottom of the conduction band 
or the top of the valence band in the insulator are defined 
as the injection barriers A„ and A p for electrons and 
holes, respectively. These barriers relate to the difference 
between the electrode work functions Ea as well as to the 
insulator interband gap energy E g : 



E A>S -A- = E A , C -A+., 



A+ 



p 



E n 



(10) 



(11) 



From now on, the minus and plus superscripts denote 
the barriers at the interfaces x = —L/2 and x = L/2, 
respectively. 

Let us characterize the insulator by the DOS functions, 
g„ t i (E) and g p j (E), describing extended states in the 
conduction and valence bands in which electron and hole 
transport takes place, respectively. Due to the large value 
of the energy gap, it is possible to introduce two sepa- 
rate electrochemical potentials; K n> i for electrons and n Pt i 
for holes. Introducing band edges means that the DOS 
function g n ^ (E) = when E < and the DOS func- 
tion g P: i (E) — when E > 0. Hence, the densities of 
electrons and holes in the extended states can be calcu- 
lated using Boltzmann statistics, thus allowing solely for 
a proper description of non-degenerate insulators, 



(x) 



[3ttV(x)] 2/3 -ecf>{x)+E b , (7) 



2m* 



where n c (x) is the electron density in the electrode, m* 
is the effective mass in the conductor and E^ is the bot- 
tom of the conduction band. Performing calculations in 
the similar manner as for the hole-injecting electrode, we 
obtain the electric field in the conductor electrode, 



F c (x) = ^ 



exp 



x - L/2 
Itf, c 



(8) 



with the unknown electric field F c (L/2) at the insula- 
tor/conductor interface again being determined by the 



n i [x) = / g n ,i [E - A„ - Koa,s + ecj) (x)J /„ (E) dE, 

(12) 



Pi ( x ) = I 9 P ,i [F + A p - Koo jS + e<j) (x)\ f p (E) dE, 

(13) 



/„ (E) = exp 



(x)-E 



kT 



(14) 
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f p (E) = exp 



E — Kp,i (x) 

kf 



(15) 



here T is the absolute temperature and k is the Boltz- 
mann constant. 

The electrochemical potentials n, n ^ and n p ^ are conse- 
quently expressed in terms of the charge-carrier densities 
rii and p h 



K n ,i (x) — feTln 



rij (x) 
N 



e<t>(x), (16) 



steady state, the total current remains constant through 
the entire system. 

Finally, from Eqs. I|21jl -(j24 |l we obtain the system of 
equations for the charge-carrier densities and the electric 
field, governed by Gauss law, in the insulator: 



kT 

n't (x) Ft (x) + n, {x) F[ (x) + — < (x) 

e 



B 

H"n,i 



R [m (x) p l (x) - n? nt ] , 



(26) 



K p,i ( x ) = —kT la 



V% (x) 
P 



- A: 



where the quantities 



ect>{x), (17) 



N = J g n (E r ) exp 
o 



EL 
kf 



dE' 



(18) 



kT 

V \ (x)Fi(x)+ Pi (x)F!(x) p'[{x) 

e 



B 



R [m (x) pi (x) - nf nt ] , 



F[ (x) = [pi (a;) - n, (x)] . 



(27) 
(28) 



and 



u 

P= j g p (E')eiq> (j£)d& 



(19) 



can be understood as the effective total densities of states 
available in the conduction and valence bands of the in- 
sulator, respectively. 

The total current density in the insulator consists of 
the electron and hole contributions and reads: 



3 = j P 0) + Jn (x) , 



(20) 



with 



jn = <Jn,i (x) Fi (x) + kT^n.iTl'i (x) , (21) 

jp = <r P ,i (x) Ft (x) - kTfip^p't (x) , (22) 

where a n ^ = e\x n ^n (x) and a Pt i = efx Pl i P (x) are the 
conductivities of electrons and holes, respectively, with 
the respective electron and hole mobilities fi ni i and fi p ^. 
When accounting for the electron-hole recombination in 
the insulator, the steady-state continuity equations for 
holes and electrons read 



dj n jx) 
dx 



eBR [m (x) pi (x) - nf at ] , 



(23) 



= -eBR [m (x) pi (x) - nf nt ] , (24) 

where < R < 1 is the recombination parameter, B is 
the Langevin recombination coefficient, 



e 

B = (fJ<n ,i + fJ*p,i) i 



(25) 



and the intrinsic charge density n- m t is given by nf nt = 
iVPexp (— Eg/kT) 18 ^ 3 . We note one more time that, in 



C. Self-consistency and boundary conditions at the 
contacts 

Assuming no dipole layer or surface charge at the in- 
terfaces one has to require continuity of the electrical 
displacement and of the electrochemical potential 40 , 



eF (x) = continuous, 
k (x) = continuous. 



(29) 
(30) 



This continuity may be provided self-consistently match- 
ing the expressions of the electrical displacements and 
electrochemical potentials. In particular, one can write 



e s F s 



e l F i 



eiFi 



L 



(31) 



Matching the electrochemical potentials we obtain four 
nonlinear boundary conditions: 



In 



In 



rij {-L/2) 
N 



p l (-L/2) 



kT kT 



-F 



3_ 

CT. 



p 
kT 



el 



TF,i 



kT 



e» \ 2 a s 



-0, 
(32) 



= 0, 
(33) 



In 



In 



(L/2) 



N 



ft (L/2) 



A+ _ el T F,c 
kT kT 



P 



A+ 
p 

kT 



elTF.c 

kT 



e M( L - 

e, \2 



-F, 1 L 



3_ 

a. 



3_ 



= o, 

(34) 

= 0, 
(35) 
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which depend on parameters of both the insulator and 
the electrodes. Notice, that from Eqs. (f3"2"|) and (|3"5|) the 
following relations are obtained: 



n % (-L/2) Pi (-L/2) = n % (L/2) Pi (L/2) 

E, 



NP exp 



J 9 

kT 



(36) 



As the last boundary condition, we can use the value of 
the current density taken at one of the interfaces: 



j=j p (-L/2)+j n (-L/2). 



(37) 



III. PHYSICAL AND NUMERICAL ANALYSIS 



in Eqs. (|34|) - ([35|) . Hence, the boundary conditions do not 
depend directly on j in most practical cases. Under this 
approximation, equations (f3"3"|) - ([3"4")) can be reduced to 



Pi 1 - ¥ 



= Pexp 
= iVexp 



p 
kT 

A+ 

n 

' kT 



el 



TF,: 



kT 

d T F,c e 
kT e. 



■F, 



, (38) 
(39) 



As a consequence, the dependence of Fi(=pL/2) as well 
as pi(—L/2) or nt{L/2) on the current is only due to 
Eq. (|3"T|) . Apparently, the injection barriers effectively 
vary with oc Fj(=fL/2) similarly to the case of unipolar 
carrier injectio n 23 ! 24 . This change in injection barriers 
leads to the definition of effective injection barriers, 



The full set of nonlinear differential equations (|26|) - 
with the nonlinear boundary conditions (|32[) - (135p 
and (|37p has to be solved numerically. 

As an example for an insulator, an organic semicon- 
ductor can be considered, exhibiting many typical char- 
acteristics of insulators such as relatively large band gaps 
up to 3 eV and, hence, the absence of intrinsic charge car- 
riers. In simple organic light-emitting diodes (OLEDs) a 
thin layer of a organic semiconductor is contacted with a 
low work-function metal and a high work-function trans- 
parent conducting oxide. Here, we consider indium tin 
oxide (ITO) as the hole-injecting electrode and Ca as the 
electron-injecting contact. ITO is typically employed as 
anode in OLEDs, since it provides a decent conductivity 
and a sufficient high work function (up to 5 eV) to allow 
for efficient hole injection while being transparent in the 
visible region of the light spectrum to ensure light out- 
coupling. On the other hand, Ca with a work function of 
about 2.9 eV delivers efficiently electrons and, hence, is 
used as cathode in OLEDs. From now on, it is assumed 
that the material specific quantities of the organic semi- 
conductor and the electrodes adopt the typical values 
given in Table |TJ 

The injection barriers A Pi „ are given by the energetic 
difference between the highest occupied molecular or- 
bital (HOMO) or the lowest unoccupied molecular or- 
bital (LUMO) in the organic semiconductor and of 
the respective electrodes. Change of the barrier heights 
at the interfaces and of the gap energy in the organic 
semiconductor, leaving the electrodes unchanged, can be 
understood as considering different organic semiconduc- 
tors. Here, for simplicity and without loss of generality, 
we assume that the gap energy in organic constituent is 
fixed. Therefore, changing the value of the barrier height 
at one interface results in equally shifted energies of the 
HOMO and LUMO levels which entails the change of the 
barrier at another interface. 

Notice, that specifying the material parameters for the 
insulator and the electrodes leads to substantial conse- 
quences for the boundary conditions. The current den- 
sity j is multiplied by the small factors iTF.s/kT^sUoc^ 
in Eqs. (f3"2")) - (133p and by the small factor ^tf.c/^TMc^oo.c 



A eff = A 



V 



" ~2 ' ' 1 



A eff 



A+ 



' F - ' ^ I Itf, 



= A+ 



( \ I 



(40) 
(41) 
(42) 
(43) 



The changes in the effective injection barriers correspond 
to the amount of energy a charge carrier gains (or loses) 
in the electric field prevailing in the electrodes. This 
change, as is seen from Eqs. (J4TJ) and (|43|) . is not neces- 
sarily negative but can be positive, since a space-charge 
region at the interface might serve as a potential barrier 
itself. The electric field, however, cannot be arbitrary 
large. Assuming Boltzmann statistics for both injected 
electrons and holes requires that the electrochemical po- 
tentials ([in]) and (fT7|) do not approach the respective 
band edges so that inequalities n, <C N ans pi <C P hold. 
Considering boundary conditions (f3"8")) - (f3"9")) . this imposes 
a requirement on the electric field at the interfaces so that 
both effective barriers (|4"0|) and (|4"2"|) remain positive. To 
keep the validity of Boltzmann statistics the above re- 
quirements are controlled during the further numerical 
calculations. 



Spatial distributions of charge carriers and 
electric field 



First, we calculate the spatial distributions of charge 
carriers and the electric field in the organic layer with- 
out accounting for recombination (R = 0). The barrier 
height for hole injection was varied from to 0.57 eV. 
For the parameters from Table HI using equations (fTU]) 
- ([TT]) and due to the fixed gap energy of 2.4 eV a hole 
injection barrier of 0.57 eV at the anode corresponds to 
a barrier-free electron injection from the opposite elec- 
trode. In Fig. [3 the distributions of holes and electrons 
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TABLE I: Typical material parameters for electrodes and an organic semiconductor (Refs. 21 41 42 43 44). The parameters are 
deduced assuming T — 300 K and m e is the free electron mass. 



ITO 


Organic semiconductor 


Ca 


Itf.s "co,s e s Ms m s Koo.s Ea,s 
(A) (cm" 3 ) f^j (me) (eV) (eV) 


N,P U fXp,i /J,n,i E g 

<™-*> (£) (£) <« v > 


ItF,c ^<oo,c ftc ^oo,c Ea,c 

(A) (cm" 3 ) (SE!) (eV) (eV) 


8.6 10 20 9.3 3 0.35 0.22 5 4.7 


10 21 3 10~ 4 10" 6 2.4 


0.8 2.6 ■ 10 22 66.7 4.68 2.87 



in the organic layer are depicted for two different cur- 
rent densities of j = 2 and 100 mA/cm 2 . The elec- 
tron and hole density profiles change slightly with the 
current but depend strongly on the respective injection 
barrier heights. When one injection barrier is small, a 
high charge carrier density appears near the respective 
electrode decreasing monotonously through the organic 
layer and dropping quickly near its ejecting electrode. In 
this case the injection barrier at the other contact is large 
and thus, the density of the respective charge carriers is 
relatively small. Due to the absence of recombination 
all injected charge carriers have to traverse the entire or- 
ganic layer. The large amount of space charge emerging 
at the low injection-barrier electrode prevent the ejection 
of the traversed charge carriers injected at the high bar- 
rier contact. This gives rise to a local maxima of their 
concentration near their ejecting electrode. The effect is 
specific as long as no hole-electron recombination occurs. 

The distribution of the electric field across the organic 
layer is shown in Fig. [3] for the same parameter values 
as in Fig. [2] The electric field strongly varies with the 
change of the injection barrier. When the heights of 
the injection barriers at the contacts differ significantly, 
one kind of the charge-carrier dominates. As a result, 
the electric field changes monotonously between the elec- 
trodes, just like in the single- interface case (see, for exam- 
ple, Ref . |24[) . In the limit of barrier- free injection the field 
profile beyond the virtual electrodeSi is well described in 
the frame of the space-charge approximation. Once the 
heights of electron and hole injection barriers at the cath- 
ode and the anode are comparable, nearly equal hole and 
electron densities prevail in the organic layer. Hence, the 
charge density is small and the electric field appears to 
be almost constant with a small maximum in the middle 
of the organic layer. 

Charge recombination gives rise to severe changes in 
the charge carrier density and field distributions in the or- 
ganic layer. First, let us consider the case of comparable 
injection barriers, i.e. A~ — 0.3 eV and A+ = 0.27 eV. 
The dependences of charge-carrier densities on coordi- 
nate x are shown in Figs. SJa and b for different val- 
ues of the recombination parameter R. The annihila- 
tion of traversing holes with traversing electrons leads to 
a strong spatial decrease in the charge-carrier densities. 
The change occurs already for a small increase of the re- 



combination parameter up to R = 0.1. A subsequent in- 
crease of the recombination parameter up to 1 influences 
the picture insignificantly. In Fig. |4jc the space distri- 
bution of the recombination rate RriiPi is shown. To il- 
lustrate the position of the recombination zone, which is 
proportional to the product riiPi, the distribution of riiPi 
normalized by its maximum value is shown in Fig. |4jd. 
For small R the recombination zone extends almost over 
the entire organic layer, but with increase of R it moves 
close to the electron-injecting electrode. This is because 
of the assumed large difference (two orders of the magni- 
tude) in the mobilities of the holes and electrons in the 
organic material. The injected holes can traverse nearly 
the entire organic layer before they recombine with the 
electrons. The influence of recombination on the respec- 
tive electric field distribution is displayed in Fig.O Even 
small values of recombination parameter (up to R = 0.1) 
change the electric field significantly, moving its maxi- 
mum to the region of the maximum of the recombination 
rate. At low R the charge density in the center of the de- 
vice is small giving rise to a nearly constant electric field. 
At the contacts however, where one charge carrier specie 
dominates space-charge effects emerge. Upon increased 
R nearly the entire organic layer is dominated by holes 
and only close to the cathode a substantial electron den- 
sity emerges. Hence, space-charge effects are important 
in the entire device. 

In the case of strongly different injection barriers their 
influence becomes more pronounced. This is illustrated 
in Figs. [6] and where the recombination zone and the 
distribution of electric field are shown for the barrier 
heights of A" = 0.4 eV and A~ = 0.2 eV, respectively. 

For A" = 0.4 eV the electric field distribution is dom- 
inated by the buildup of negative space charge. Due to 
the high hole mobility injected holes traverse the organic 
layer and recombine close to the cathode where the max- 
imum electron density prevails. This holds for all R and 
thus, the recombination maximum is situated close to 
the cathode. For A" = 0.2 eV, when the injection of 
holes is more efficient than the injection of electrons, the 
field distribution is dominated by the buildup of positive 
space-charge. For very small R the maximum of recom- 
bination zone is situated near the hole- injecting electrode 
since the injected electrons can traverse the whole organic 
semiconductor without recombining and their density is 
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still locally increased near the hole-injecting electrode. 
However, with increased R the recombination probabil- 
ity increases and the few injected electrons recombine 
instantaneously after injection since their recombination 
rate exceeds their transport rate. Hence, the maximum of 
the riiPi product moves to the cathode with increased R. 

B. Current-voltage characteristics 

Next, we calculate the current-voltage characteristics 
of the system under consideration. Knowledge about the 
distribution of the electric field gives access to the voltage 
drop V across the system for a given current density j 
and hence, to its IV characteristics. 

The voltage drop at the device is defined as 

/+oo 
F d (x)dx-V bi , (44) 
-OO 



where 

( F s (x)-j/a s , x < -L/2, 

F d {x) = { Fi (x) , -L/2 <x< L/2, (45) 
[ F c {x) - j/a c , x > L/2, 

and — Vbi is the voltage drop in the case of j = 0, i.e. 
the built-in potential. Using the system of equations 
(|2l) ]l -(|25 j) and boundary conditions (|g2" ]) -([55 ]) . the built-in 
potential can be obtained analytically. This calculation 
is presented in the Appendix. Similarly to the case of 
unipolar charge-carrier transport, the built-in potential 
is given by the difference of the electrode's work func- 
tions, 

eV bi =E AtC -E AtS . (46) 
Integrating the field Fd, we obtain 



Fd (x) dx = I 



TF,: 



-Fi 



J_ 

0". 



+ 1 



TF.i 



V 



—Fi [ L 
2 



L/2 



Fi (x) dx. 



(47) 



-L/2 



Using the boundary conditions, the voltage may be 
rewritten in the shorter form 



V 



kT . 



■In 



p l {-L/2) 



Pi {L/2) 



L/2 



Fi (x) dx. 



(48) 



-L/2 



In Fig. [5] the resulting IV characteristics are presented 
for the case of absence of recombination when the bar- 
rier height Ap changes from to 0.57 eV. At voltages 
V < — Vu there is no Ohmic-like j ^ V behavior which 
has been obtained in the cases of the unipolar transport 23 
or of the bipolar charge injection from Ohmic contacts^. 
In the calculated range of voltages (and currents) all IV 
characteristics exhibit only a tendency to this behav- 
ior. Hence, it can be supposed that the Ohmic region 
lies at even lower voltages. Near the built-in voltage, 
— Vu, a relatively wide region starts where, for the cases 
of barrier- free injection of carriers (holes or electrons), 
the SCLC behavior j ~ V 2 of IV characteristics occurs. 
This behavior is similar to that described in the diffusion- 
free analysis in Ref . and stems from the dominance of 
electrons or holes, respectively. When nonzero injection 
barriers exist for both kinds of charge carriers, the IV 
dependences differ from the SCLC behavior but all cor- 
responding curves lie between the curves with A~ = 
and Ap = 0.57 eV. For high voltages at all A p the IV 
dependences tend to a transition in the j ~ V 3 regime 
which is intrinsic for systems with bipolar transport in 
the recombination-free case^i. Yet, in the framework of 
our model this region can not be reached at some barrier 



heights because the used Boltzmann statistics of carriers 
in the organic constituent is violated here. 

When recombination is accounted for, the IV charac- 
teristics change substantially. This influence is demon- 
strated in Fig. [§] on the example of characteristics for 
the barrier heights of A~ = 0.3 and 0.4 eV. One can see 
that the most significant effect occurs in the high- voltage 
part of the characteristics where the transition from the 
j ~ V 3 regime to the SCLC-like j ~ V 2 behavior oc- 
curs even for the smallest value of R = 0.001. At the 
same time, the remaining part of the IV characteristics 
virtually is not changed by the recombination. However, 
it should be noted that with increase of R the current 
density slightly increases in the low voltage part of the 
IV characteristics whereas in the region of the SCLC-like 
behavior a slight decrease of j is visible. 

The significant change of the IV characteristics at high 
voltages can be related to the influence of the recom- 
bination on the modification of the injection barriers. 
The corresponding dependences of the effective barriers 
A g |j on the applied voltage V are presented in Fig. [TO] 
In general, the effective barriers decrease with the volt- 
age tending to become equal at vanishing barrier for the 
recombination-free case. However, recombination weak- 
ens this decrease so that nonzero effective barriers are 
obtained at much higher voltages. By that the difference 

between A~rr and A + „. remains considerable and even 
_efi_ eft 

increases significantly for high magnitudes of the recom- 
bination rate. Paradoxically, the prevailing injection bar- 
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riers provide a significant extension of the voltage range 
where the SCLC-like behavior of the IV characteristic 
occurs. 



C. Recombination current and efficiency 

integrating the product eBR [n< (x)pi (x) — nf Dt \ [i.e. 
the right-hand side of Eq. (p?3")) ] over the thickness of the 
organic layer, we obtain a quantity having the dimension 
of a current density and usually called the recombina- 
tion current density, j r (Ref. 0). It stands for the total 
number of recombination events in the volume of organic 
constituent per unit square of the interface, per unit time. 
This quantity determines the recombination efficiency of 
the device, r\ r = j r /j (Ref. 0). The calculated depen- 
dences of the recombination current density on the ap- 
plied voltage are given in Fig. QT] for the barrier heights 
Ap = 0.3 and 0.4 eV. The total recombination current 
increases with increased recombination parameter R and 
thus a higher luminance of the OLED is expected. It 
is clearly seen from Fig. [11] that already for the small- 
est calculated values of the recombination parameter R 
a substantial recombination current is observed while in 
the range from R = 0.1 to 1 the increase in j r rather 
weak. 

In general, the recombination currents follow the IV 
characteristics with an increased total recombination rate 
for an increased current density. However, the depen- 
dence of the recombination efficiency rj r on the applied 
voltage unveil important details. In Fig.[T2]the recombi- 
nation efficiency is depicted using the same parameters 
as for the system in Fig. 111! The recombination efficiency 
shows features which have not been discussed in previous 
works 2 i 29 i 33 i 34 i 35 ] 36 where different boundary conditions 
at the injecting interfaces have been assumed. In gen- 
eral, a maximum in r\ r (V) is found close to —Vu no mat- 
ter which injection barriers have been considered. This 
may be explained by the detailed analysis of the carrier 
densities and partial currents as follows: a fast increase 
of the electric current j occurs at voltages V < — Vu 
(Ref. [23J) corresponding to the Shockley diode equation 
for unipolar injection^. Assuming a virtually indepen- 
dent unipolar injection for electrons and holes from the 
two electrodes one can expect the faster increase of the 
recombination current j r quadratic in carrier densities 
which results in the peak close to —Vu- At higher bias 
r) r (V) depends strongly on the injection barrier heights 
and the charge carrier mobilities. For rather balanced 
but high injection barriers (A~ = 0.3 eV) there is a large 
range of voltage where the efficiency is strongly attenu- 
ated independently of the recombination parameter. As- 
suming that the built up of space-charge is small due 
to the substantial injection barriers, the electric field is 
constant in the device and the injection barrier lowering 
scales linearly with the applied voltage as can be seen 
from Fig. [TUla in between 3 and 30 V. However, for low 
bias the injection barrier lowering is negligible so that the 



density of electrons and holes stays roughly constant for 
a wide voltage range. Therefore the recombination cur- 
rent does not change if the bias increases. At the same 
time, since the total current is proportional to the elec- 
tric field a drop of the efficiency upon the bias increase 
is expected. For an initiation of barrier lowering, the 
charge carrier densities of electrons and holes grow expo- 
nentially with the barrier lowering and thus an increase 
in the recombination current sets in. The efficiency sat- 
urates and increases again. For a further increase in bias 
the reduction of the injection barrier heights is weakened 
since space charges emerge. The respective injection bar- 
riers tend to diverge leading to a higher injection barrier 
for holes than for electrons balancing the charge carrier 
densities in the insulator. Since balanced charge-carrier 
densities result in an optimized recombination rate the 
self-balancing effect leads to an increase in efficiency. 

Introducing asymmetric injection barriers the 
efficiency-voltage characteristic changes dramatically 
For a higher value of the hole injection barrier, 
Ap = 0.4 eV, and thus a low electron injection barrier 
the attenuation of the efficiency due to high injection 
barriers is missing. The most pronounced change in the 
efficiency curve is the appearance of the peak denoted 
by D in Fig. [T^b for R = 1. Increasing the voltage 
coming from —Vu leads to a buildup of space-charge 
for the charge carriers with the lower barrier, in this 
case electrons. Strong injection of electrons contributes 
to the electric current j but does not add much to the 
recombination current because the minority carriers 
remain in the barrier-limited regime. This provides the 
fall of the efficiency down to the point C. However, since 
the electron mobility is low, the negative space charge 
does not result in a strong current. On the other hand, 
the Langevin recombination parameter B is large since 
it is dominated by the larger hole mobility. Therefore, 
the efficiency is not fully attenuated. In the area C 
to D where the electrons are already in the space-charge 
regime with approximately constant and low barrier, 
the barrier for holes decreases (see Fig. [TUlb) providing 
exponential increase of the minority carriers. This 
does not contribute much to the electric current j but 
strongly promotes the recombination current j r . Further 
increase of the voltage depresses the injection barrier at 
the hole-injecting side leading to an substantial amount 
of holes accumulated at the anode. Now the above 
mentioned self-balancing effect sets in and the efficiency 
tends to increase again. Apparently, the charge carrier 
mobilities are important for the efficiency peak D. 
Assuming comparable charge carrier mobilities for 
electrons and holes, a strong reduction of the efficiency 
can be observed below 30 V. 

It should be also noted that the change of the mobil- 
ity ratio (for example, by decrease of the hole mobility) 
not only varies quantitatively the values of the current 
and the efficiency but also acts strongly on the position 
and the shape of the recombination zone. This is seen 
from Fig. 1131 where the spatial distributions of this zone, 
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calculated in several points of the r\ r (V) dependences of 
Fig. [121 b for the cases [i v ^ — n n: i and fi Pt i — 100 n n ,i, are 
presented. With increased applied voltage, the recombi- 
nation zone extends through the whole organic layer as 
long as equal mobilities fi Pj i — fi n ^ are assumed, whereas 
in the case of the large mobility difference there is a pro- 
nounced peak of the recombination zone moving to the 
electron-injecting electrode since the injected holes are 
faster than the electrons. 

Thus, additionally to the case of fully balanced injec- 
tion and transport properties, which is often difficult to 
achieve in single-layer devices^, it may be possible to 
maximize the recombination efficiency in an OLED with 
imbalanced injection and transport properties compen- 
sating the small amount of the minority carriers by their 
relatively high mobility with respect to the mobility of 
the dominating carriers. 



D. Comparison with experiment 

To evaluate the presented model it is compared to ex- 
perimental data obtained from a diode consisting of a 
single poly(phenylene vinylene) (PPV) layer of thickness 
L = 100 nm sandwiched between ITO and Ca electrodes. 
The characteristic energies of the LUMO and HOMO 
levels of the employed PPV are 2.8 eV and 5 eV, respec- 
tively and the hole mobility fj, p j equals 5-10 -7 cm 2 /(V s) 
(Refs. 6.32). The ratio of the electron and hole mobilities 
is assumed to be [i n ,i/[i<p,i = 0.01. It is known, that the 
work function of ITO is sensitive to the cleaning proce- 
dure and thus, it may be varied from 4.7 eV to 5 eV (up 
to the energy of the HOMO level). Therefore, we are al- 
lowed to consider the injection barrier A~ as one of the 
fitting parameters. The other barrier, A+, is supposed 
to be unchanged and equals 0.07 eV. 

In Fig. [14] the measured IV characteristic of the 
ITO/PPV/Ca structure is shown as well as the IV de- 
pendences calculated for A~ = 0.1 eV (the best fit- 
ting value) and for different values of fitting parameters 
R and N(= P). On the one hand, the calculation satis- 
factory reproduces the magnitude and the general form 
of IV characteristic using reasonable parameters. On the 
other hand, however, the exact shape of the IV curve can 
not be fully approximated by a unique set of parame- 
ters. This is due to the fact that the proposed model 
misses still some important features of insulators, in par- 
ticular of organic semiconductors. First of all, realistic 
DOS shapes for the LUMO and HOMO levels as well as 
trap levels in the insulator layer have to be accounted 
for. This leads effectively to a field and charge carrier 
density dependence of the mobilities which may increase 
the calculated currents at larger voltages. It should be 
also noted, that due to this the ratio of the hole and 
electron mobilities varies with the change of the applied 
voltage^. 

Comparing the luminance efficiency in Fig. [T5] with the 
calculated recombination efficiency a good reproduction 



of the voltage dependence can be obtained above 3 V. 
Here, the calculated recombination efficiency was nor- 
malized to the measured luminance efficiency ignoring 
e.g. the voltage dependence of the extraction efficiency. 
There is, however, a strong discrepance of the measured 
and fitted efficiency curves at lower voltages. Again, this 
difference may be due to the absence of trap levels in our 
simple model. At low voltages, strong trapping could 
arise close to the electrode shifting the recombination 
zone to the injecting contact giving rise to a strong re- 
duction of the luminance efficiency due to quenching ef- 
fects&M 

Concentrating on the effect of injection on the device 
performance one has to consider also a possible forma- 
tion of the dipole layers at the electrode/organic interface 
which may have crucial consequences for the prevailing 
charge-carrier injection barriers. Dipole layers may be a 
result of self-organization or a consequence of the depo- 
sition of some functional layer at the interface^^. The 
phenomenon of a dipole-layer formation cannot be incor- 
porated self-consistently in the developed phenomenolog- 
ical approach, however, an effective change in the injec- 
tion barrier heights can be assumed. Notice, that in this 
case the magnitudes of injection barriers A~ and A+ are 
decoupled and become independent. 

If compared to other models of bipolar transport and 
recombination in OLEDs the model presented here de- 
scribes satisfactorily a crossover from the high-voltage, 
space-charge limited regime to the low-voltage, barrier- 
dominated regime. Even the regime around the built- 
in potential which is problematic in theories using other 
boundary condition a 21 ' 32 i 48 is well described. Comparing 
with very comprehensive discrete one-dimensional simu- 
lations by Tutis et al. 37 our approach presents a simpler 
continuous description using a strongly reduced number 
of fitting parameters. 



IV. CONCLUSIONS 

In the calculation of the charge-carrier transport 
through insulators, the fundamental question about 
boundary conditions generally arises when charge- 
carrier injecting interfaces are involved. In this pa- 
per, a simple one-dimensional model describing the 
bipolar charge-carrier transport across a semiconduc- 
tor/insulator/conductor structure was presented, where 
the problem of injection was for the first time solved 
self-consistently. Essentially, continuity of the electric 
displacement and the electrochemical potential was as- 
sumed, thus matching them at the injecting interfaces 
of the electrodes with the insulator. Additional bound- 
ary conditions were defined far into the electrode materi- 
als, where the influence of the involved materials on each 
other is negligible so that they can be regarded as inde- 
pendent. Considering the Poisson equation, the electric- 
field distributions and the current-voltage characteristics 
were derived using as an example the organic semicon- 
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ductor sandwiched between ITO and Ca electrodes. The 
influence of different injection conditions and recombi- 
nation rates on the current transport properties of such 
a structure was analyzed. It was found that the injec- 
tion barriers determine the dominating carrier type and 
strongly influence the distributions of the charge and of 
the induced field in the organic layer. The recombina- 
tion influences on these quantities most strongly when 
the heights of injection barriers for holes and electrons 
are comparable. It was shown that the injection barriers 
and recombination together with the mobility ratio de- 
fine the position and shape of the recombination zone in 
the organic layer as well as drastically affect the current- 
voltage characteristics of the considered diode structure. 
We have also found that in the case of significantly dif- 
ferent injection barriers a high recombination efficiency 
of a diode may be still achieved in the case of a relatively 
high mobility of the minority carriers. Notice, that the 
controlling of the carrier mobility in organic devices is 
experimentally possible^. 

Finally, we can conclude that the presented self- 
consistent model of charge injection reveals unambigu- 
ously that injection barriers and recombination rate are 
the governing factors controlling the transport properties 
of semiconductor/insulator/conductor structures. On 
the other hand, this model cannot yet provide a com- 
plete description of the particular case of OLEDs because 
some specific features of these systems are still missing. 
Most important of them are the realistic DOS shapes 
including energetically distributed trap state s 14 ' 30 ' 31 ' 32 
or concentration- and field-dependent carrier mobilities 50 
characteristic of organic semiconductors. The latter fea- 
ture is followed by the exponential-type dependence of 
the injection current on the device thicknes s 36 ' 51 which 
can hardly be expected within our consideration. Fur- 
ther extension of the model concerning the inclusion of 



the above-mentioned features is necessary for better ap- 
plication to experimental data and is now in preparation. 
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APPENDIX: CALCULATION OF THE BUILT-IN 
POTENTIAL IN THE CASE OF BIPOLAR 
TRANSPORT 

The built-in potential defines as 

/+oo 
Fd(x)\ j=0 dx, (A.l) 
-oo 

with Fd (x) from Eq. ([45)) . In the case of j = both 
current densities j p (holes) and j n (electrons) are equal 
to zero. From Eqs. [|2T]) - (p2]) we get 



kT 

— raj (x) + m (x) Fi (x) = 0, 

e 

kT 

— v\ (x) - Pi (x) Fi (x) = 0, 



or, combining these expressions, 

„ . „ kT d , 
Fi (x = — — In 
Ze ax 



Pi (g) 

7li (x) 



dx 



In \pi (x) rii (x)] = 0. 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



Integrating in Eq. (|A.1[) with Fi (x) from (|A.4[) . we find 



J 



kT 

V bl = Itf,sF s {-L/2) + l TF , c F c {L/2) + — 

Ze 



In 



Pi (L/2) 
i (L/2) 



n 



In 



Pi (-L/2) 



n t (-L/2) 



r 



(A.6) 



Using the boundary conditions ([32|) -([35 | with j = 0, we 
directly obtain 



Vu = -(K-K. 



(A.7) 



Accounting for Eq. (fl0|) . it follows immediately that 
eVbi=E AtC -E A)S . (A. 



It should be also noted, that from Eq. (|A.5|) follows 
that the product pi (x) rii (x) takes the constant value 
everywhere in the organic layer. To provide zero current 
in the absence of applied voltage and be consistent with 
Eqs. (|26ll - (f27|) . we choose this constant equal to nf nt . 
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FIG. 1: Schematic band diagram of the considered semiconductor/insulator/conductor structure. 
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FIG. 2: Spatial distributions of the charge-carrier densities, and m, in units of no = EieokT/e 2 l^ F s for a system without 
recombination for different barrier heights A~ and a constant current density: [(a) and (b)] j = 2 mA/cm 2 and [(c) and (d)] 
j = 100 mA/cm 2 . 
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FIG. 3: Spatial distributions of the electric field Fi in units of Fo = kT /elrF,s for a system without recombination for different 
barrier heights A~ and a constant current density: (a) j = 2 mA/cm 2 and (b) j = 100 mA/cm 2 . 
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FIG. 4: Spatial distributions of (a) the hole density pi, (b) the electron density rii, (c) the recombination rate, and (d) the 
recombination zone ~ riiPi for a system with the barrier height A~ = 0.3 eV and the current density j = 2 mA/cm 2 for 
different values of the recombination parameter R. The carrier densities are measured in units of no = ueokT/e 1 T f,s- The 
product riiPi in the part (d) is normalized on its maximum value. 



~° 0.08 
ll" 0.04 



R = 


R = 0.1 


R = 0.001 


R = 1 


" R = 0.01 




* **" ' 


' N 




\ 


j = 2 mA/cm 


■ ■•■ ■ 






A = 0.3 eV; A + = 0.27 eV 

p n 
i . i . i . i 





FIG. 5: Spatial distributions of the electric field Ft in units of Fo = kT/elrF,s for a system with the barrier height A~ = 0.3 eV 
and the current density j = 2 mA/cm 2 for different values of the recombination parameter R. 
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FIG. 6: Spatial distributions of (a) the recombination zone ~ riiPi and (b) the electric field Fi for a system with the barrier 
height Ap — 0.4 eV and the current density j = 2 mA/cm 2 for different values of recombination parameter R. The electric 
field is measured in units of Fq — kT/elrF,s- The product iiiPi is normalized on its maximum value. 
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FIG. 7: The same as in Fig. [6] but for the injection barrier height A p = 0.2 eV. 




FIG. 8: IV characteristics of a system without recombination for barrier heights A p = 0, 0.2, 0.3, 0.4, and 0.57 eV. 
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FIG. 11: Dependence of the recombination current density on the voltage V for a system with barrier heights A p = 0.3 eV 
(solid curves) and 0.4 eV (dotted curves) for different values of the recombination parameter R. 
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FIG. 12: Dependence of the recombination efficiency r) r on the voltage V for a system with barrier heights (a) A~ = 0.3 eV and 
(b) 0.4 eV for different values of the recombination parameter: R = 0.001 (squares), R = 0.1 (circles) and R = 1 (triangles). 
The ij r (V) dependence for A~ = 0.4 eV, equal carrier mobilities fi P; i = fi n ,i and R — 1 is also shown by dashed curve. 
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FIG. 13: Spatial distributions ol the recombination zone ~ riiPi at different voltage for a system with the barrier height 
Ap = 0.4 eV, the recombination parameter R — 1 and the hole mobilities fj, Pi i — fj, n< i (upper four figures) and [i v ^ = 100 fi n ,i 
(bottom four figures) for the points of the r/ r (V) curves marked in Fig. 1121 b by capital letters from A to D. The carrier densities 
are measured in units of no = eie^kT / e 2 l^ F . 
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FIG. 14: Measured IV characteristic of ITO-PPV-Ca structure (solid curve) and its approximation with different values of the 
model parameters (dashed curves with symbols). 




FIG. 15: Measured voltage dependence of the luminance efficiency of ITO-PPV-Ca structure (solid curve) and its fitting with 
the model parameters (dashed curve). The last dependence is shown in arbitrary units. 



